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ABSTRACT 

In a recent paper we presented the calculation of the 0{as) virtual 
corrections to 6 — s and of those bremsstrahlung terms which are 
needed to cancel the infrared divergences. In the present paper we work 
out the remaining 0{as) bremsstrahlung corrections to 6 — >■ sl'^l~ , which 
do not suffer from infrared and collinear singularities. These new contri- 
butions turn out to be small numerically. In addition, we also investigate 
the impact of the definition of ttIc on the numerical results. 
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1 Introduction 



The inclusive rare decay B ^ Xg l'^t~ has not been observed so far, but is expected to be 
measured at the operating B factories after a few years of data taking. The measurement 
of its various kinematical distributions, combined with improved data on i? — 7, will 
imply tight constraints on the extensions of the standard model and perhaps even reveal 
some new physics. 

The main problem of the theoretical description of i? — > is due to the long-distance 

contributions from cc resonant states. When the invariant mass ^/s of the lepton pair is 
close to the mass of a resonance, only model dependent predictions for these long distance 
contributions are available today. It is therefore unclear whether the theoretical uncertainty 
can be reduced to less than ±20% when integrating over these domains |^. 

However, when restricting to a region below the resonances, the long distance effects are 
under control. The corrections to the pure perturbative picture can be analyzed within the 
heavy quark effective theory (HQET). In particular, all available studies indicate that for 
the region 0.05 <s = s/ml < 0.25 the non-perturbative effects are below 10% 0, ^, |], ^, 
^, |]. Consequently, the differential decay rate for B ^ Xg can be precisely predicted 
in this region using renormalization group improved perturbation theory. It was pointed 
out in the literature that the differential decay rate and the forward-backward asymmetry 
are particularly sensitive to new physics in this kinematical window |T0|, |ri|, p!3|] . 



The next-to-leading logarithmic (NLL) result for B ^ Xg £^£ suffers from a relatively 



large (±16%) dependence on the matching scale [14, 15|. The NNLL corrections to the 
Wilson coefficients remove the matching scale dependence to a large extent , but leave 
a ±13%- dependence on the renormalization scale /ij,, which is of 0{mh). In order to further 
improve the theoretical prediction, we have recently calculated the 0{as) virtual two-loop 
corrections to the matrix elements {s i~^i~\Oi\b) {i = 1, 2) as well as the virtual 0{as) one- 
loop corrections to O7,..., Oio [0. As some of these corrections suffer from infrared and 
collinear singularities, we have added those bremsstrahlung corrections needed to cancel 
these singularities. This improvement reduced the renormalization scale dependence by a 
factor of 2. 

In the present paper we complete the calculation of the bremsstrahlung corrections as- 
sociated with the operators Oi, O2, 07,...,0io, i.e., we add those bremsstrahlung terms 
which are purely finite and have therefore been omitted in Ref. [0]. We anticipate that the 
additional terms have a small impact on the phenomenology of 6 — > s£^i~. 

The paper is organized as follows: In Sec. |^, we briefly specify the theoretical framework, 
before, in Sec. |^, we discuss the organization of the calculation of the finite bremsstrahlung 
corrections and review the structure of the virtual corrections and singular bremsstrahlung 
contributions, calculated in Ref. [Q. The finite bremsstrahlung corrections are worked 
out in Sec. § and Sec. In Sec. |, finally, we investigate the numerical impact of the 
new corrections on the invariant mass spectrum of the lepton pair. We also illustrate the 



dependence of our results on the definition of the charm quark mass. 



2 Effective Hamiltonian 

The appropriate tool for studies on weak B mesons decays is the effective Hamiltonian 
technique. The effective Hamiltonian is derived from the standard model by integrating 
out the t quark, the Zq and the W boson. For the decay channels b ^ s {i = /i, e) it 
reads 

Heff = — -y=- v^i Vtb y ^ Cj Oj , 

where Oi are dimension six operators and Cj denote the corresponding Wilson coefficients. 
The operators we choose as in [ p!6[] : 



Oi = 




O2 = 


(sl7mCl)(cl7^&l), 


03 = 




O4 = 




05 = 




Oe = 




Or = 


Js 


Os = 




O9 = 




Oio = 





The subscripts L and R refer to left- and right-handed fermion fields. We work in the 
approximation where the combination {Vus^uh) of Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements is neglected and the CKM structure factorizes. 

In the following it is convenient to define the operators O7,..., Oio according to 

0, = ^0,, (J = 7,...,10), (1) 

with the corresponding coefficients 

C, = —C,, (j = 7,...,10). (2) 



3 Organization of the calculation and previous results 

In this section we comment on the organization of the calculation of the virtual and brems- 
strahlung corrections to the process b ^ s and repeat the results obtained in Ref . [|l| . 

The one- loop diagrams in Fig. ^ associated with the four-quark operators Oi,..., Og, lead 
to contributions which are proportional to the tree level matrix elements of the operators 
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b Ci Oi s b Cj Oi s 




Figure 1: Diagram (a) can be absorbed by replacing the Wilson coefficients Cj and Cg 
through and C^°^, respectively. 7* denotes an off-shell photon which subsequently 

decays into a pair. Similarly, diagram (b) is absorbed through the replacement 

Cg — >■ C™°*^. g denotes an on-shell gluon. The index i runs from 1 to 6. See text for 
details. 








Figure 2: Diagrams which are automatically taken into account when calculating correc- 
tions to ''""'^^Oy, Ct'^'"'^d, and ct'^'"'^d,. 



O7, Os and O9. JTherefore, they can be absorbed by appropriately modifying the Wilson 
coefficients C7, and Cg. The modified coefficients we write as 



^mod 

8 



(3) 



: Ag + Tg h{z, S) + [/g s) + Wg /l(0, s) 



The auxiliary quantities Aj, Tg, t/g and Wq are linear combinations of the Wilson coefficients 
Ci{n). Their explicit form is relegated to the appendix. The one-loop function h{z,s) is 
given by ||16| 



h{z, s) 



4, , , 8 16 z 

--Inf^) H \ 

9 ^ ^ 27 9 s 



2 + — 




Az - s 



2 arctan 



iz—s ' 



In ( ^+^^ 



I TT, 



s <Az 



s> Az 



(4) 



It is obvious that the modification of the Wilson coefficients automatically accounts also 
for the diagrams in Fig. ^ when calculating the corresponding corrections to the matrix 
elements 



I a 



(O,mod) 



7,8,9), 



3 




Figure 3: The two-loop virtual diagrams induced by Oi and O2 that cannot be absorbed 
into the 07,8,9 contributions by weighing them with the modified Wilson coefficients. The 
circle-crosses denote the possible locations where the virtual photon is emitted. The curly 
lines represent gluons. 




Figure 4: The only two bremsstrahlung diagrams induced by Oi and O2 that cannot 
be absorbed into the 07,8,9 contributions by weighing them with the modified Wilson 
coefficients. 



where are the leading order terms of the modified Wilson coefficients, i.e., 

p;(0,mod) _ .(0) 
L/y — /I7 , 

^(O,mod)^^(0)^ (5) 
^(O,mod) ^ ^(0) ^ y(0) ^^^^ ^ ^(0) ^^^^ ^ ^(0) ^^Q^ 

For the explicit expressions of the quantities Af'\ Tg°\ Ug'^^ and Wg^^ we refer to the 
appendix. 

Notice that the virtual and bremsstrahlung corrections of the four-quark operators with 
topologies shown in Figs. ^ and ^, however, have to be calculated explicitly. As the Wilson 
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a) b) c) d) e) 

Figure 5: One-loop virtual 0{as) corrections induced by c!f'^°'^^Oj, C^'^°'^^Os, 6*9°'™°'^^ O9 
and C^Oio- The circle-crosses denote the possible locations for emission of a virtual 
photon. 




a) b) c) 



Figure 6: The 0{as) bremsstrahlung diagrams induced by O7, O9, Oio and Og. Weighing 
the contributions of O7, Og and Og with the corresponding modified Wilson coefficients 
accounts for the bremsstrahlung diagrams depicted in Fig. ^ (b)-(e). The crosses and 
circle-crosses denote the possible locations for emission of a bremsstrahlung gluon and a 
virtual photon, respectively. 



coefficients Ci and C2 are much larger than C3, 
topologies only for Oi and O2 insertions. 



,C6 we retain the contributions of these 



, we systematically calculated the virtual corrections to the matrix 



In the previous work 

elements of C[^^Oi, C^'02, shown in Fig. ^, as well as to those of C'^^'^°"^'Oj {j = 7, 9) 

and C[^Oio (cf Fig. Furthermore, we also took into account the corrections to the 
Wilson coefficients calculated in Refs. 



16, ff^ 



We found that the matrix elements of the operators O7, Og and Oiq [cf Fig. ^(a)-(c)] 
suffer from infrared and collinear singularities. Consequently, on decay width level the 
interferences {Oj,Ok) {j,k = 7,9,10) are singular, too. We therefore included the gluon 
bremsstrahlung corrections associated with {Oj,Ok) {j,k = 7,9,10) in order to get an 
infrared finite result for the decay width [cf Fig. ^(a) and (b)]. 

Taking into account the virtual and bremsstrahlung contributions discussed so far, we 
obtain the result presented in Ref. ||l||: 



dr{b^Xsi+t 



ds 



'l + 2s) 



CefF 
9 



+ 



48 7r3 

2^ 



:i-s) 



+ 4(l + 2/s) cf + 12 Re ( Cf Cg"'^* 



(6) 



where the effective Wilson coefficients Cy , Cg and Cf^ are given by 0] 



=( l + ) Aj 



TT 



{cr^Fr\s)+Cr'FPis)+APFPis)), (7) 



Cf =fl + ^fM u;g{s) \ {Ag + Tg /i(m2, §) + Ug h{l, s) + Wg h{0, s)) 



4:71 



(cf (s) + cf Ff (5) + Af Ff (S) 



C^it= ( 1 + ^^9(5) ) ^10. (9) 



TT 



The quantities Cf , Cf , ^7, /if, Ag, Aio, Tg, f/g and Wg are Wilson coefficients or 
linear combinations thereof. We give their analytical expressions and numerical values 
in the appendix. The one-loop function h{ml,s) is given in Eq. (^), while the two-loop 
functions f'^2^'^\ accounting for the diagrams in Fig. and the one- loop functions F^^'^^\ 
corresponding to the diagrams ^(d) and (e), are given in Ref. The functions oj-j and cjg, 
finally, include both virtual and bremsstrahlung corrections associated with O7, Og and 
Oio- For details on their construction we again refer to [|l|. 

When calculating the decay width (^, we retain only terms linear in as (and thus in cuy, cug) 
in the expressions for jC^'^^p, |Cg^p and |C^oP- ■'■^ interference term Re [Cf^Cf^* ) too. 



we keep only linear contributions in as- By construction one has to make the replacements 
cug ujjg and —>■ ujjQ in this term. 

The functions uj, ug and turg read 

18 + s, 2 s(2 -2s - s^) , 1 16 -lis -1752 

ln(l - s) - - ^ ^ ln(s) 



3 2 + s ' ' 3(l-s)2(2 + s) 18 (2 + s)(l-s) 



4 2 2 5 + 4 s 

^9(5) = -g Li(s) - - ln(l - s) ln(s) - -tt^ - ^^g. ln(l - s) 



(10) 



2s(l + g)(l-2g) 5 + 9g-6g^ 

3(1- s)2(l + 2 s) ^ 6 (1 - s)(l + 2 s) ' ^ ^ 
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9 s ^ ' 9 ^ ^ 18 1 - s ^ ' 



Summary 

The bremsstrahlung corrections associated with the interferences 



(Cf '"^"^^O,, d^'^^'-'^O,) , (j, k = 7,9, 10), 



are aheady included in formula @. The remaining bremsstrahlung corrections, which are 
infrared finite, we derive in Sec. ^ and Sec. ^ In Sec. ^ we discuss the contributions of the 
interferences 

which we call to be of type A. There is no contribution from k = 10 because of the Dirac 
structures of the involved operators. Sec. |^ is devoted to the interferences 

(Cf )0„ Cf O,) , it, J = 1, 2) and (cf 0„ Cj^'^^^^d,) , = 1, 2; A; = 7, 8, 9, 10) . 

Accordingly, we call these the type B terms. Again, the contributions for = 10 vanish 
due to the Dirac structures of the operators involved. 



4 Finite bremsstrahlung contributions of type A 

The bremsstrahlung contributions taken into account by introducing the functions uJi{s) 
cancel the infrared divergences associated with the virtual corrections. All other brems- 
strahlung terms are finite. This allows us to perform their calculation directly in c? = 4 
dimensions. 

The bremsstrahlung contributions from O7 — Og and Os — Og interference terms as well as 
the Os — Os term oppose no difficulties. The sum of these three parts can be written as 

^pBremSjA ^pBrems ^pBrems ^pBrems 

78 _j_ 89 _j_ 88 



ds ds ds ds 

) ) 48^^^ I ^"""^ + ^''^ + )■ ^^^^ 

The coefficients Cij are given by 

2 



C78 — Oi? ■ Og , Cgg — Oi? ■ Og Og , Cgg — Up ■ Og , [i^j 



^(o,cfr) 



7 



while the quantities Tij read 



= — <^25 - 27r^-27s + 3s^-s=^ + 12 (s + s^) ln(s) 
9s * 



TT 



+ 6 I — — arctan 



2 - 1 ,S + 



(2 - s)V^V4-s 
12 ( (1 - s)VI "\/4 - s - 2 arctan 



24 Re I^Li 

V4 - s 



2-s 



X 



arctan 



4- s 



arctan 



\/4 - .s 
2-s 



7"88 



- Stt^ + (1 - s) (77 - s - 45^) - 24Li(l - s) 



27 s 

+ 3|10-4s-9s^ + 81n 





) 


l-s 





ln(s) +48Re|^Li 



3-s .(l-s)V4^ 



-6 



20s + 10s2-3s3 



— 8 TT + 8 arctan 



4- s 



X 



arctan 



4- s 



— arctan 



2-s 



T89= 3<|5(4-s)-3-41n(s)(l 



- 8 Rc ^Li 



s . V§ -\/4 — s 
2+^ 2 



Li 



-2 + s(4-s) .(2-s)v^V4^ 
;7 + I' —i:: 



+ 4 I a/ + 2 arctan 



arctan 



2-s 

— arctan 



X 



4- s 



2-s 
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5 Finite bremsstrahlung contributions of type B 



In this section we consider the bremsstrahlung contributions from Oi and O2 and interfer- 
ence terms with O7, Og, O9 and Oiq. As mentioned before, interferences with Oiq vanish 
due to the Dirac structures of the operators. 

The bremsstrahlung contributions discussed in this section all involve the matrix elements 
associated with the two diagrams depicted in Fig. ^. Their sum, J^p, is given by 



T 6 9s Qu 

16 Tl'^ 



E{a, P, r) A«5 + E{a, P, q) Ai^ - E{P, r, q)—Ai2z 

q-r 

-E{a, r, q)^ Azae - E{/3, r, q)^ Ki^j 
q-r q-r 



l\. (18) 



where q and r denote the momenta of the virtual photon and of the gluon, respectively. 
The index a will be contracted with the photon propagator, whereas (3 is contracted with 
the polarization vector e^(r) of the gluon. Jq,/3 and Kik are obtained from and /S.ik 
im], respectively, by setting = and dropping terms proportional to rp. The matrix 
E{a,j3,r) is defined as 

E{a,p,r) = ^(7a7/3/ - flpla)- (19) 
Due to Ward identities, the quantities Ai^ are not independent of one another. Namely, 

q°'Jaf3 = and r^Jap = 
imply that Ai^ and Aig can be expressed as 

- - - 

A25 = Ai23 + — Ai27 ; Aie = Ai26 • (20) 

q-r 

As in addition Az26 = — A«23, the bremsstrahlung matrix elements depend on Az23 and 
Az27, only. In c? = 4 dimensions we find 

C 

Az27 = 8(g-r) / dxdy ^^^'^^ , (21) 



Ai23 = 8 (q-r) / dx dy 
Jo 



where 

C = ml — 2xy{l — y){q-r) — q^ y {1 — y) — i6. 

In the rest frame of the b quark and for fixed s = q'^/ml, the phase space integrals which 
one encounters in the calculation of dT^^^'^'^'^ / ds can be reduced to a two-dimensional 
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integral over Er — Er/rrii, and Eg — Eg/mb, where Er and Eg are the energy of the gluon 
and the s quark, respectively. In the following it is useful to introduce the integration 
variable w — 1 — 2Es instead of Eg. The integration limits are then given by 



Er e 



w — s w — s 



2w 



and w e [s, 1] 



For fixed values of s, the quantities Ai23 and Ai27 depend only on the scalar product q-r, 
which, in the rest frame of the b quark, is given by {w — s)ml/2. The integration over 
Er turns out to be of rational kind and can be performed analytically. The remaining 
integration over w, however, is more complicated and is done numerically. The result can 
be written as 



^pBrems,B 

ds 



Utt/ Utt/ 



^b,pole \ Ws ^tfeP 

487r3 



X 



± 

J dw (Cii + Ci2 + C22) T22 + 2 Re [ (Ciy + C27) T27 + (Ci8 + C28) T28 + (cig + C29) T29] 



(22) 



The quantities Tij, expressed in terms of Ai23 and Ai27, read 
8 (w-s)(l-wy 



T22 = 



27 



X 



3w^ + 2r(2 + w) -sw(5-2w) 



2 8^(2+ w) +sw(l + 2w) 



|Ai27| +4s 



w {1 — w) — s {2 + w) 



|Ai23| + 

Re [Az23Az;,] I (23) 



8 1 

T27 = o~ 

3 s w 



1^1 — w) (4 — s w + w^) + s w (4 + s — w) ln( 



Ai 



23 



4 (1 - w) + s w (4 + s - w) ln(w) Ai27 ^ (24) 



^28 = 7T~ 



9 s w (w — s) 



X 



[w - sf{2s-w){l-w) 



23 



2s(^/;-s)^(l-w) Ai27 



+ sw 



(1 + 2 s - 2 w)Ai23 - 2 (1 + s - w)Ai27 



■In 



1^1 + s — w){w'^ + s (1 — w)) 



(25) 



4 1 
6 w 



2 s{l — w){s + w) + 4:swln{w) 



Ai 



23 



2 s(l - 'u;)(s + w) + ^(3 s + w) ln{w) A^ar ^ (26) 
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The coefficients q-,- in Eq. ( p2[ ) include the dependence on the Wilson coefficients and the 
color factors. Explicitly, they read 



Cll = Cr 



(0) 



^ _r r^{o)/=;(o,efr)* 



Ci2 = Cro ■ 2 Re 



C22 — Cp 



a 



(0) 



2 



7 ; 
(0)p;(0,eff)* 



C18 — Ct2 ■ C*! 
^ _ ^(0) P;(0,eff)* 



C27 — Lyp -0207 , 

C28 — tvF ■ tv2 *-"8 ' 

^ ^{o)p;{o,cfr)* 

C29 — L^F ■ '-"2 '-"9 5 



(27) 



where the color factors Cp, and C^-j arise from the following color structures: 



^T'^T'^ = Cf1, Cf 



2Nr 



a,b,c 



8iV| 



and 



a,b 



Am ■ 



Finally, we list the explicit formulas for Az23 and Az27 expressed as a function of s and the 
integration variable w. We obtain 



Ai23 = -2 + 



— s 



(28) 



AZ27 = 2 



(29) 



where the functions Gk{t) {k > —1) are defined through the integral 



Gk{t)= I dxx" \n[l-tx{l-x) -i6], Gi{t) = -Go{t). 
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Explicitly, the functions and G'o(^) read 

'2 7r arctan(^y^j - ^ - 2 arctan^J^y^ j , t<A 

-2z7rln(^^±^) -^ + 2 1n2(^^^*±fOj, t>4 



(30) 



TT,/^ -2-2 arctanW^ j , t<4 



-t7rJ'-^-2 + 2Ji^ ln( ^+f^ j , t>A 



(31) 



6 Numerical results 



First, we investigate the numerical impact of the finite bremsstrahlung corrections [see 
Eqs. (|T^ and (P2|)] on the dilepton invariant mass spectrum. Following common practice, 
we consider the ratio 



quark ^ 



dT{b^si+£~) 



T{b — > XcCZ/e) 



ds 



(32) 



in which the factor ml^^^y^ drops out. The explicit expression for the semi-leptonic decay 
width r(6 — >• XcCUe) reads 



Tib^X,eu,)= ^^^^3 iKfcl ^7|-2 



c,pole \ ^/^c 



6, pole 



1 2 I ' 



where (/(z) = 1 — 8 z + 8 — z'^ — 12 z"^ ln{z) is the phase space factor, and 



ir(^) 



2as{mb) f{z) 



37r 5((z) 



25 239 25 
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z + —z']+z \Yi{z) 20 + 90 2 - - 2^ + — 2 



17 64 



17 



+ z' \n\z) (36 + z') + {l-z') ^-~^-z + - z' j ln(l - ^) 
- 4 (1 + 30 + z^) \n{z) ln(l - 2) - (1 + 16 + z*) (6 Li(2) - vr^) 



32z^/^{l + z) 



TT 



4 Li(v^) + 4 Li(- v^) - 2 ln(z) In 







Vi + 





(33) 



(34) 



incorporates the next-to-leading QCD correction to the semi-leptonic decay |jT8[. The 
function f{z) has been calculated analytically in Ref. |19|. It reads 



(35) 
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0.05 
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Figure 7: The new contribution Ai?quark(S) due to finite bremsstrahlung corrections for 
/i = 2.5 GeV (uppermost curve), = 5 GeV (middle curve) and /x = 10 GeV (lowest 
curve) and rric/mb = 0.29. 

We stress that the function f{z) refers to on-shell renormalization of the charm quark 
mass. 

In Fig. I^we consider the contribution Ai?quark('S), which is due to the finite bremsstrahlung 
corrections in Eqs. ([I3|) and (p2D, for three values of the renormalization scale (/i=2.5, 5 and 
10 GeV) and for fixed value mc/m^ = 0.29. The values of all the other input parameters are 
as in Ref . [|^ . In Fig. ^ we combine the new corrections with the previous results. The solid 
lines show the ratio -Rquark(s), including the new corrections, for the values /i = 10 GeV 
(uppermost curve), /i = 5 GeV (middle curve) and fi = 2.5 GeV (lowest curve) and for fixed 
value rric/mb = 0.29. The dashed lines represent the corresponding results without the new 
corrections. We find that for s = 0.05 the new corrections increase the ratio -Rquark(S) by 
~ 3%, while for larger values of s their impact is even smaller. When including the finite 
bremsstrahlung corrections we obtain 



for the integrated quantity -Rquark- The error is obtained by varying /i between 2.5 GeV 
and 10 GeV. For comparison, the corresponding result without the finite bremsstrahlung 
correction is i?quark(s) = (1.25 ± 0.08(/i)) x 10"^ §]. 

Among the errors on -Rquark(s) due to uncertainties in the input parameters, the one related 



0.25 



'quark 




,-5 



0.05 



13 




0.25 



Figure 8: The solid curves show the ratio -Rquark(s) including the finite bremsstrahlung 
corrections while the dash-dotted curves show the corresponding results without the new 
corrections. The uppermost curves (solid and dash-dotted) correspond to /i = 10 GeV, the 
middle curves to fi = 5 GeV and the lowest curves to yU = 2.5 GeV. rric/mb = 0.29. 



to the charm quark mass is by far the largest. We therefore only comment on this error. In 
principle, the uncertainties induced by the charm quark mass have two sources. First, it is 
unclear whether ttLc in the virtual- and bremsstrahlung corrections should be interpreted 
as the pole mass or the MS mass (at an appropriate scale). Second, the question arises 
what the numerical value of rric is, once a choice concerning the definition of rric has been 
made. 

To illustrate these problems more clearly, it is useful to first consider the process B Xs'y. 
There, the one-loop matrix elements of Oi and O2 vanish, implying that the charm quark 
mass dependence only enters at 0{as)- Formally, one can interpret rric in these 0{as) 
expressions to be the pole mass or the MS mass because the difference is of higher order 
in ag. Nevertheless, it has been argued in the literature [2^] that the choice m^^(/i) 
with /i G [rric, f^b] seems more reasonable than m^"^^ (which was used in all the previous 
analysis) due to the fact that the largest charm quark mass dependence comes from the real 
part of the two-loop matrix elements of Oi and O2, where the charm quarks are usually 
off-shell, with a momentum scale set by ml"^^ (or some seizable fraction of it). It was 
shown in Ref. that the definition of the charm quark mass leads to a relatively large 



uncertainty in the branching ratio: Changing rric/mh in V[B — > Xs^) from 0.29 ± 0.02 to 
0.22±0.04, i.e., from mP°^'^/m{j'°''' to m^^/m^"^*^ (with /i G [nic, rab\), causes an enhancement 
of BR{B^Xs-i) by ~ 11%. 
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0.06 0.08 0.1 0.12 0.14 0.16 
s 

Figure 9: i?quark(s) for various values and definitions of rric: Tlie three bands are obtained by 
setting mP°^^/m^°^'^=0.31 (uppermost), 0.29 (middle) and 0.27 (lowest) in T(b XcCVe)- 
In the rare decay h XJ+l' we set mf^/mP°^'' = 0.18,0.22,0.26. This leads to three 
curves all within a narrow band. See text. 



In the process B ^ Xs this problem is less severe because vtlc enters already the 
one-loop diagrams (i.e., at 0(a°)) associated with 0\ and O^- As the two-loop calculation 
requires the renormalization of mc, the definition of nic has to be specified. Therefore, the 
two-loop result explicitly depends on the definition of the charm quark mass. This can be 
seen from For the pole mass definition, the results for the two-loop matrix elements of 
0\ and O2, encoded in F^'^2^^\ are given in Eqs. (54)-(56), while those corresponding to 
the MS definition are obtained by adding the terms AF{^*2mcren given in Eq. (49). 

In the following, we investigate the impact of pole- vs. MS definition of rric in the rare 
decay b ^ Xs on the ratio -Rquark("S)- In the semileptonic decay b ^ X^eUe the charm 
quark is basically on-shell. Therefore, we always use the pole mass definition for the charm 
quark mass in T{b — > XcCZ/g), which enters -Rquark('5)- 

In Fig. H we set mP°^'^/m^°^'^ equal to 0.31, 0.29 and 0.27 in the decay width T{b XcC z/g). 
In the rare decay b ^ Xs ^ on the other hand, we use the MS definition for mc, and put 
^MS^^poie _ 0.18,0.22 and 0.26 (independently of mP°'^/m^°^'^, to be on the conservative 
side). This leads, for a given value of mP°'^/m^°'^, to three curves which form a narrow 
band. The uppermost band corresponds to mP°^'^/m^°^'^ = 0.31, the middle to 0.29 and 
the lowest to 0.27. The curves with the strange behavior for s > 0.13 all belong to the 
lowest value m^^/m^°'^ = 0.18. As the result for the two-loop corrections was derived in 
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Figure 10: i?quark(s) for various values and definitions of mc- The solid curves are obtained 
by setting mP°^^/m^°^'^=0.33 (uppermost), 0.31, 0.29, 0.27 and 0.25 (lowest) in the rare- 
and the semileptonic decay. The dashed lines are obtained by taking m^^/m^°^'^ = 0.22 in 
the rare decay and mP°^'^/m^°^'^= 0.31, 0.29 and 0.27 in r{b ^ X^e z/g). See text. 



expanded form which only holds for s < 4m^/m^, the strange behavior illustrates that, 
for mc/rrih = 0.18, the result is not valid for s > 0.13. 

In Fig. |1^ the three middle solid curves are obtained by adopting the pole mass definition of 
nic, both in the rare and in the semileptonic decay. They correspond to mP°'''/m^°'^ = 0.31, 
0.29, 0.27. The dashed curves, on the other hand, are obtained when the MS definition 
with m^^/ml°^'^ = 0.22 is used in the rare decay width. One finds that for s > 0.06 the 
results for -Rquark("S) are somewhat larger when using the pole mass definition of rric in 
the rare decay. For values below s < 0.06 the situation is reversed and thus the same as 
for b Xs'~f ^0|1- Again, the strange behavior of the dashed curves indicates that, for 



rric/mb = 0.22, the expanded formulas become unreliable for values of ,§ > 0.19 . The thick 
solid lines are obtained by adopting the pole mass definition on the whole and correspond 
to mc/m(,=0.33 (upper) and 0.25 (lower). In summary, the figure shows that the quark 
mass uncertainties can effectively be estimated by working with the pole mass definition 
throughout, provided one takes the rather conservative range 0.25 < mP°i7mP°^' < 0.33. 

Finally, in Fig. |IT] we show -Rquark('5) in the full range s G [0.05,0.25] for mP°^'^/m^°^° G 
[0.25,0.33]. Note that for these values of mc/mb the expanded formulas hold just up to 
s = 0.25. 

Comparing Fig. |] with Fig. we find that the uncertainty due to rric/ is clearly larger 



16 



Q I I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 
§ 

Figure 11: i?quark(s) for mP°^7mP°^^=0.33 (uppermost), 0.31, 0.29, 0.27 and 0.25 (lowest) 
in the rare- and the semileptonic decay in the full window s e [0.05, 0.25]. 



than the leftover fi dependence. Varying rric/mfj between 0.25 and 0.33, the corresponding 
uncertainty amounts to ±15%. 

To conclude: We have calculated the finite gluon bremsstrahlung corrections of 0{as) to 
r(6 — > si~^i~), taking into account the contributions of the operators Oi, O2, O7, Oq, Og 
and Oio- We have worked out the numerical impact of the new corrections on the invariant 
mass spectrum of the lepton pair in the range s e [0.05, 0.25] and found an increase of 
about 3% for s = 0.05 and even less for larger values of s. Furthermore, we investigated 
the uncertainties of i?quark(s) due to the definition and numerical uncertainties of the 
charm quark mass. We found that these errors can be reliably estimated when working 
with the pole mass definition of rric, provided one takes the rather conservative range 
0.25 < mP°i7mP°^" < 0.33. 
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A Auxiliary quantities Ai, Tg, Ug and Wq 



The auxiliary quantities Ai, Tg, Ug and Wg, appearing in the modified Wilson coefficients 
in Eq. (|^) and in the effective Wilson coefficients in Eqs. are the following linear 

combinations of the Wilson coefficients Cj(/i) |T2|: 



A, 
As 

Ag 

Aio 

Tg 
Ug 
Wg 



Air 



20 



80 



9 



9 



^ Csifx) + Csifx) - I Cifx) + 20 C5(/i) - ^ C6(/i) 



«s(/i) 



6 



C^io(/i) 



(36) 



- C^{^^) + C2(/i) + 6 CM + 60 CM , 



7 



2 32 

-CM'^^CM-^CM 



^CM)- 

12 32 

o CM - o ^4(/i) - 8 CM - -T CM ■ 



The elements 7|g'* can be found in |]T6[, while the loop-function h{z, s) is given in Eq. (^. 

In the contributions which explicitly involve virtual or bremsstrahlung correction only the 
leading order coefficients Af\ Tlf \ u!f^ and Wg^^ enter. They are given by 



A 



(0 



A 



(0 



A 



(0 



C 



(1) 

7 



1 



(1) 



Cs 

in 

a. 



C 



cP 



(0) 



6 ' 



4 ^(0) 20 (0) 80 (0) 
~9^' "Y^' 



(0) 



47r 



6 " ' 3 



rrib 



4(0) 


— '-'10 ' 










+ 6Cf 


+ 60 Cf ^ , 




""2^=^ "3 


^ - 3 


o ^(0) 32 (0) 

^'-^S ""3" ^6 




""2^=^ "3 


) - 8 


(0) 32 (0) 

'-'5 ^'-"6 • 



(0) 



64 



^5 



64 



.(0) 



(37) 



We list the leading and next-to-leading order contributions to the quantities A^, Tg, Ug and 
Wg in Tab. § 
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2.5 GeV 


5 GeV 


10 GeV 


as 




0.267 


0.215 


0.180 






-0.697 


-0.487 


-0.326 






1.046 


1.024 


1.011 






(-0.360, 0.031) 


(-0.321, 0.019) 


(-0.287, 0.008) 






-0.164 


-0.148 


-0.134 


(4'. 




(4.241, - 0.170) 


(4.129, 0.013) 


(4.131, 0.155) 


(rf. 




(0.115, 0.278) 


(0.374, 0.251) 


(0.576, 0.231) 


(cf. 




(0.045, 0.023) 


(0.032, 0.016) 


(0.022, 0.011) 






(0.044, 0.016) 


(0.032, 0.012) 


(0.022, 0.009) 






(-4.372, 0.135) 


(-4.372, 0.135) 


(-4.372, 0.135) 



Table 1: Coefficients appearing Eqs. (0)-(P) for /x = 2.5 GeV, /x = 5 GeV and /i = 
10 GeV. For as{n) (in the MS scheme) we used the two- loop expression with five fiavors 
and as{mz) = 0.119. The entries correspond to the pole top quark mass = 174 GeV. 
The superscript (0) refers to lowest order quantities. 
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